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from  disabled  high-speed  aircraft.  These  aircraft  have  been  in  existence 
for  over  three  decades  and  there  have  been  numerous  elections  bv  nilots  of 


many  countries.  In  these  thirty  years,  there  has  also  been  remarkable 
advances  in  aircraft  and  space  technologv.  However,  progress  in  the  area 
of  injury  reduction  in  pilot  ejection  has  been  slow.  One  of  the  major 
factors  for  the  lack  of  progress  is  the  assumption  that  the  spine  is  re- 


quired to  carry  the  inertial  load  of  the  torso,  proportional  to  the  ac- 
celeration level  of  the  seat,  and  that  there  is  very  little  anvone  can 
do  to  reduce  this  load  for  a given  acceleration.  This  concept  is  quite 
evident  in  the  early  models  for  simulating  the  ejection  event.  Latham  (1) 
assumed  a lumped  parameter  base  excitation  model  in  which  the  torso  or 
mass  was  placed  above  a spring  representing  the  snine.  The  assembly  was 
accelerated  upward  by  a force  or  acceleration  pulse  applied  to  the  bottom 
of  the  spring  and  the  resulting  peak  force  in  the  spine  was  greater  than 
the  product  of  mass  and  acceleration.  It  was,  in  fact,  dependent  upon 
the  rate  of  onset  of  the  input  pulse.  Conceptually , this  model  provided 
an  over-simplification  of  the  problem  and  stifled  anv  ideas  of  load  re- 
duction borne  by  the  spine.  The  continuum  model  proposed  by  Hess  and  Lom- 
bard (2)  was  also  too  simple,  since  it  simulated  the  acceleration  of  a 
straight  rod  or  column  which  had  a uniform  weight  distribution  along  its 
length. 

Subsequent  models  by  many  investigators  perpetuated  the  axial  loading 
concept  and  reduction  of  injury  was  not  achieved  through  these  efforts. 

In  1968,  King  et  al  (1)  proposed  and  verified  the  hypothesis  that  the  spine- 
al  column  carried  an  eccentric  load,  and  that  it  was  subjected  to 
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a simultaneous  compressive  and  bending  load.  The  basis  for  this  hypothesis 
was  strain  gage  data  from  instrumented  vertebrae  of  embalmed  cadavers 
which  were  subj ected  to  +G^  impact  accelerations,  using  a vertical  accel- 
erator at  ;Wayne  State  University. 

Dr.  C.  L.  Ewing  of  NAMRL  had  been  interested  in  the  pilot  ejection 
problem  and  had  observed  that  anterior  wedge  fractures  were  a common  injury 
mode  in  pilots  who  eject.  He  proposed  a hypothesis  that  the  posterior 
structure  of  the  spine  acted  as  a motion  limiter,  preventing  spinal  hyper- 
exter.sion,  but  facilitating  spinal  flexion.  Ewing  proposed  to  verify  this 
hypothesis  by  putting  the  spine  into  moderate  extension  to  position  the 
vertebrae  in  a less  favorable  position  for  anterior  wedging.  It  was  rea- 
soned that  posterior  compression  of  the  vertebral  column  is  limited  tv  the 
articular  facets  while  anterior  compression  is  unlimited. 

In  view  of  the  observed  bending  effects  and  the  Ewing  hypothesis, 

i 

a research  project  was  conceived  to  verify  the  hypothesis  and  to  study 
the  precise  mechanisms  of  vertebral  injury  due  to  +0^  acceleration.  The 
proposed  sponsor  was  the  Office  of  Naval  Research.  The  specific  aims  of 
the  study  were: 

1.  To  verify  that  spinal  hyperextension  will  reduce  strain  along  the 
anterior  aspect  of  vertebrae  in  the  thoraco-lumbar  spine. 

2.  To  determine  the  extent  and  manner  of  hyperextension  which  will 
produce  the  maximal  reduction  in  strain. 

3.  To  verify  that  there  is  a statistically  significant  increase  in 
fracture  g-level  sustained  by  cadaveric  subjects  as  a result  of 
spinal  hyperextension. 

Upon  sucessful  attainment  of  these  aims,  the  following  objectives  were 
pursued:  , 


.. 
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A.  To  quantify  the  existence  of  a facet  lead  which  constituted  a 
second  load  path  along  the  vertebral  column. 

5.  To  identify  the  role  of  abdominal  pressure  as  a possible  third 
load  path. 

6.  To  determine  the  effect  of  the  spinal  musculature  on  vertebra] 
column  response  due  to  +G^  acceleration. 

7.  To  develop  a two-dimensional  mathematical  model  of  the  spine 
simulating  its  response  to  +C^_  acceleration. 

8.  To  conduct  cadaveric  experiments  which  will  provide  data  for  the 
validation  of  this  model. 

9.  To  propose  this  model  for  the  assessment  of  injury  probability 
in  ejection  seat  design. 


II.  THE  EWING  HYPOTHESIS 


Several  hypotheses  have  been  advanced  to  explain  these  vertebral  frac- 
tures. Most  of  the  proposed  explanations  concern  simply  peak  acceleration 
values,  and  peak  rates  of  onset  of  acceleration,  and  vie"  the  vertebral 
column  as  a single  structure  or  series  of  single  structures  having  a 
physical  failure  point,  which,  when  exceeded,  causes  structural  failure 
of  one  or  more  vertebrae. 

It  is  hypothesized  by  others,  however,  that  such  structural  failure 
limits  can  be  increased  by  altering  body  positioning.  In  this  view, 
opening  the  angle  that  the  longitudinal  axis  of  the  torso  forms  with  the 
vertical  axis  of  the  unrestrained  pelvis,  by  means  of  increasing  the  seat 
angle,  will  prevent  fracture,  by  arranging  the  impact  vector  of  the  election 
thruster  so  that  it  does  not  occur  normal  to  the  superior  or  inferior  sur- 
face of  the  thoraco-lumbar  vertebral  bodies. 

Ewing  et  al  (A)  proposed  that  one  of  the  major  causes  of  ejection 


vertebral  fracture  is  the  dynamic  reaction  of  the  vertebral  column  under  +G 

z 

(eyeballs  down)  impact  acceleration  in  the  presence  of  improper  restraint; 
that  is,  there  are  certain  movements  of  the  individual  vertebral  bodies 
under  -K;  impact  acceleration  that  cause  the  characteristic  ejection 
vertebral  fracture.  If  this  theory  is  correct  and  if  these  motions  are 
prevented,  the  fractures  would  therefore  be  expected  to  occur  only  at 
markedly  higher  levels. 

This  may  be  explained  by  considering  the  vertebral  column  as  a series 
of  spring-mass  systems,  with  the  intervertebral  disks  serving  as  springs, 
the  vertebral  bodies  (and  body  segments  that  they  support)  as  individual 
masses,  and  the  anterior  and  posterior  interspinous  ligaments  as  spring  lim- 
iters (on  tension  only) • 

Figure  1 demonstrates  an  idealized  vertebral  segment  consisting  of 
three  vertebrae  acting  as  masses  and  two  intervertebral  disks  acting  as 
springs.  If  this  segment  acted  as  a simple  spring-mass  system,  an  impact 
acceleration  in  the  vector  normal  to  the  superior  or  inferior  surface  of 
the  vertebral  bodies  and  passing  through  the  center  of  mass  should  result 
in  plateau  compression  fracture(s),  as  noted  in  Figure  2,  if  the  compression 
failure  limit  of  the  superior  surface  of  the  vertebral  body  is  exceeded. 
(Acceleration  vector  direction  is  shown  in  the  figures  by  an  arrow) . Such 
fractures  are  rarely  seen  in  ejection  vertebral  fractures,  however. 

If  an  acceleration  is  applied  with  the  vector  either  before  or  behind 
the  center  of  mass,  anterior  or  posterior  wedge  compression  fractures, 
respectively,  should  result  if  the  failure  limit  of  the  superior  surface 
of  the  vertebral  body  is  exceeded  due  to  anterior  or  posterior  bending 
of  the  multiple  spring-mass  system,  as  illustrated  in  Figure  3.  however, 
recent  examination  of  eighty  ejection  vertebral  fracture  cases  by  Ewing  (5), 
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showed  that  almost  all  were  anterior  compression  fractures,  as  demonstrated 
in  Figure  4;  only  one  posterior  compression  fracture  was  found,*  One 
possible  explanation  is  that  the  posterior  compression  of  the  system  is 
subject  to  some  spring-limiting  mechanism,  whereas,  anterior  compression  is 
not.  Examination  of  the  anatomy  reveals  this  to  be  the  case. 

The  posterior  compression  limiter  (or  spring  limiter)  is  seen  to  be 
the  articular  facets  of  the  vertebrae,  held  together  by  ligaments,  as 
demonstrated  in  Figure  5,  that  serve  as  a posterior  hinge  for  adjacent 
vertebrae.  As  can  be  seen  in  Figure  6,  this  hinge  allows  anterior  vertebral 
lips  to  touch  but  prevents  any  contact  of  the  posterior  vertebral  lips.  Thus 
an  additional  spring-limiter  system  other  than  the  anterior  and  posterior 
interspinous  ligaments  is  acting. 

An  hypothesis  is  therefore  presented;  posterior  compression  of  the 
vertebral  column  in  the  thoraco- lumbar  area  is  limited  bv  the  articular 
facets  of  the  vertebrae  while  anterior  compression  is  not  limited.  If  this 
hypothesis  is  true,  a means  of  preventing  anterior  vertebral  compression 
fracture  is  suggested. 

If  the  vertebral  column  could  be  forcibly  restrained  during  application 
of  impact  acceleration  in  a position  of  relative  hyperextension  as  in  Figure 
7,  anterior  compression  would  be  limited.  Therefore,  to  cause  a compression 
fracture,  the  spinous  processes  connecting  the  articular  facets  to  the 
vertebral  body  would  have  to  be  torn.  The  forces  required  to  cause 
fractures  of  the  spinous  processes  with  restraint  in  a position  of  hvperex- 

* Autopsies  performed  on  ejection  fatalities  rarely  include  vertebral  body 
examinations.  The  conclusion  was  based  upon  examination  of  the  x-ray  re^* 
port  in  each  instance  of  those  confirmed  by  x-ray  as  having  suffered  a vert- 
ebral fracture  and  survived.  Those  fatally  injured  at  time  of  ejection  were 
not  included. 
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tension  as  in  Figure  b,  would  be  increased  quite  considerably  over  those 
required  to  cause  anterior  vertebral  fracture  while  restrained  in  the  con- 
ventional erect  manner.  Thus  the  vertebral  fracture  threshold  limit  would 
be  markedly  increased. 

III.  THE  HYPEREXTENSION  TESTS 

1.  The  Hyperextension  Device 

The  hyperextension  device  consisted  of  a wooden  block,  15  cm  x 10  cm  x 
57  cm  , that  was  fastened  to  the  seat  back.  The  height  of  the  block  from 
the  seat  pan  was  adjustable.  The  long  axis  of  the  block  was  placed  horizon- 
tally against  the  seat  back  so  that  the  block's  smallest  dimension  was  be- 
tween LI  of  the  cadaver  and  the  seat  back. 

Placement  of  the  blocks  was  determined  by  experiments  carried  out  to 
determine  the  optimal  location,  size  and  shape.  Based  on  those  findings, 
the  centerline  of  eack  block  was  placed  opposite  the  body  of  LI  in  almost 
all  cases.  Due  to  the  difficulties  inherent  in  performing  x-rays  on  a 
metal  seat  positioned  over  an  elevator  shaft,  a wooden  mock-up  seat  was 
used  to  determine  the  location  of  LI.  While  LI  was  known  to  be  instrumented 
and  this  could  be  detected  easily  on  x-ray,  the  posterior  exterior  anatomi- 
cal localization  of  LI  was  necessary  in  order  to  place  the  hyperextension 
device.  This  was  accomplished  by  using  radio-opaque  markers  on  the  wooden 
chair,  in  the  midsagittal  plane.  The  x-ray  showed  the  position  of  the 
external  markers  and  che  position  of  Ll  relative  to  them.  The  cadaver  was 
then  placed  in  the  ejection  seat,  the  geometry  duplicated  by  means  of 
identical  markers,  and  the  centerline  of  the  block  placed  next  to  the  center- 
line  of  Ll. 

2.  The  Vertical  Accelerator  Experiments 

All  experiments  were  performed  on  a vertical  accelerator  housed  in  an 
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eight-story  elevator  shaft  of  the  School  of  Medicine  at  Wayne  State  i’niver- 
sity.  The  sled  has  an  aircraft  ejection  seat  mounted  directly  above  the  pis- 
ton. The  angle  between  ceat  pan  and  seat  back  is  90  degrees.  Accelerator 
line  of  thrust  is  parallel  to  the  seat  back.  The  cadaver  was  positioned 
in  the  seat  prior  to  firing.  The  acceleration  pulse  was  induced  with  the 
sled  at  rest  at  the  bottom  of  the  shaft;  stroke  length  is  2.44  m.  At 
completion  of  the  input  acceleration  pulse,  braking  was  initiated  and 
deceleration  completed  within  about  12  m.  The  acceleration  pulse  was 
approximately  trapezoidal  in  shape,  the  rate  of  onset  and  the  magnitude 
of  the  plateau  being  variable.  Details  of  the  accelerator  have  been 
described  by  Patrick  (6). 

Foil-type  strain  gages,  0.125  in.  in  length,  we/e  bonded  to  the 
cleaned  surfaces  using  Eastman  910  adhesive  and  catalyst.  Strain  gages 
were  installed  on  the  anterior  surface  of  the  vertebra]  body  in  the  mid- 
sagittal  plane,  and  lateral  gages  were  mounted  in  bilaterally  symmetrical 
pairs  coplanar  in  the  coronal  plane  of  T12,  L2,  and  L4  on  two  of  the  cadavers. 
One  cadaver  also  had  an  anterior  gage  on  LI  and  two  posterior  gages  on  L2. 

On  the  other  8 cadavers,  the  only  gages  installed  werp  anterior  ones  on 
each  vertebra  from  T11-L4.  On  the  last  2 cadavers,  onlv  an  anterior  gage 
on  L4  was  installed,  in  order  to  preserve  the  integrity  of  the  anterior 
• spinous  ligament,  as  a control.  All  gages  were  applied  with  their  sensi- 

tive axes  parallel  to  the  vertical  axis  of  the  vertebral  body,  and  all  leads 
were  presoldered  tc  the  gage  terminals.  An  attempt  was  made  to  install 
these  gages  near  the  neutral  axis  of  the  vertebra  so  that  their  output 
indicated  predominantly  axfal  compression. 

The  gages  were  tested,  and  a minimum  of  500  Mr  of  resistance  to  ground 
was  obtained  to  ensure  long-lasting  and  noise  tree  outputs.  After  testing 
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for  proper  operation,  the  installed  gages  were  coated  with  Cagekote  it 3 to 
insure  insulation  from  body  or  embalming  fluids.  The  leads  from  the  gages 
were  attached  by  sutures  to  an  adjacent  body  part  to  prevent  accidental 
removal  during  exposure  to  experimental  conditions. 

After  the  cadaver  was  positioned  in  the  accelerator  chair,  the  leads 
from  the  lateral  strain  gages  were  connected  to  form  diagonally  opposite 
arms  of  a four-arm  Wheatstone  bridge;  the  other  two  arms  were  121-ft  high- 
stability  precision  resistors.  This  configuration  resulted  in  summing  of 
the  output  of  the  left  and  right  gages  to  eliminate  the  effects  of  lateral 
bending  of  the  vertebral  column.  The  anterior  gage  formed  a two-arm 
bridge  with  a 121-fi  resistor,  with  the  other  two  arms  being  provided  by  the 
bridge  balance  unit. 

The  following  gage  indentification  symbols  were  used:  the  prefix  A de- 
noted a gage  on  the  anterior  surface  of  T12.  The  prefix  D denoted  gages 
mounted  on  the  lateral  surface  of  the  vertebra.  Strain  values  of  DL2  re- 
presented the  average  output  of  the  two  gages  on  the  lateral  surface  of  L2 
on  either  side  of  AL2.  Posterior  gages  were  denoted  by  the  prefix  DD. 

The  restraint  system  consisted  of  an  automotive  lap  belt  under  a 
regular  aircraft  lap  belt  and  shoulder  harness,  with  leg  straps.  The  wrists 
were  tied  together  and  loosely  anchored  to  the  seat  base  without  tension 
by  means  of  a single  rope  run  through  an  eye  bolt  in  the  seat  base  to 
prevent  flailing  of  extremities.  A rope  was  also  looped  around  the  aircraft 
lap  belt  and  through  the  eye  bolt  to  serve  as  an  inverted  V.  This  step  was 
taken  so  that  no  load  would  be  transmitted  to  the  arm  rest,  and  to  prevent 
vertical  displacement  of  the  lap  belt  during  pretensioning  of  the  shoulder 
harness,  and  during  the  experimental  acceleration.  The  lap  belt  was  always 
snugly  tightened,  and  the  shoulder  harness  was  then  preloaded  to  a 89  N 
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tension  before  each  run  except  for  those  runs  in  the  flexed  mode  (vide  infra). 

Although  the  head  was  unrestrained,  its  initial  position  was  kept  approxi- 
mately vertical  bv  means  of  masking  tape  that  broke  once  the  head  started 
rotating  but  did  not  have  sufficient  breaking  strength  to  serve  as  a motion 
limiter  under  acceleration. 

3.  Fracture  Levels  and  Spinal  Modes 

Instrumented  cadavers  restrained  in  the  seat  were  subjected  to  +0 

z 

impact  accelerations  in  3 g-4  g (peak)  increments,  keeping  the  rate  of  onset 
constant  at  roughly  the  rate  of  onset  of  the  standard  eiection  seat  (approx. 

300  g/sec-500  g/sec).  The  input  acceleration  pulse  was  trapezoidal,  with 
peak  plateau  amplitudes  varying  from  5.3  g to  24.5  g.  Duration  of  the  input 
pulse  varied  from  approximately  140  m/sec  to  300  m/sec  and  was  dependent  upon 
the  peak  acceleration,  since  total  stroke  length  was  2. .44  m.  After  each  impact 
exposure,  the  cadaver  was  x-raved  again,  and  if  no  fracture  was  noted,  the 
next  incremental  exposure  was  given.  In  some  cases,  the  x-ravs  were  devel- 
oped after  the  completion  of  several  runs.  If  fracture  occurred  in  one  of 
these  runs,  the  data  from  all  postfracture  runs  were  discarded.  The  end 
point  for  each  cadaver  was  fracture  of  a vertehral  bodv  demonstrable  by 

x-rav. 

' 

At  each  peak  acceleration  level,  the  cadaver  was  tested  in  one  or  more 
of  the  following  spinal  modes:  1)  the  entire  thoraco-lumbar  segment  was 
forced  into  moderate  hvper-extension  and  maintained  in  this  position  bv  the 
restraint  system  and  hvperextei  sion  device — the  extended  mode;  2)  the  verte- 
bral column  was  allowed  to  assume  the  normal  erect  configuration  of  a 
seated  cadaver,  with  a Ian  belt  and  shoulder  harness  restraint  system — the 
erect  mode;  and  3)  the  cadaver  was  restrained  in  a seated  position  bv  a lap 


belt  only.  In  the  absence  of  a shoulder  harness,  hvperflexion  of  the 
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vertebral  segment  was  permitted — the  flexed  mode.  On  some  cadavers  anterior 
strain  gages  were  used  on  all  vertebral  bodies  from  Til  through  L 4,  but  on 
controls  only  L4  was  instrumented  to  determine  the  effect  of  instrumenting 
the  vertebrae  on  fracture  level. 

| 

The  effect  of  evisceration  on  fracture  level  and  measured  strain  was 
tested  on  the  two  cadavers  on  which  only  L4  was  instrumented  and  which  were 
not  eviscerated.  Retention  of  the  abdominal  viscera  can  onlv  affect  the 
fracture  level  adverselv  since  some  of  this  weight  must  be  borne  by  each 
spinal  segment.  To  determine  the  effect  of  the  experimental  acceleration 
as  a cause  of  fracture,  precise  fracture  timing  relative  to  onset  of  sled 
acceleration  by  examination  of  vertebral  body  strain  curves  was  used.  From 
postrun  x-ray  studies  alone  it  was  not  possible  to  determine  whether  a 
detected  fracture  was  due  to  the  experimental  acceleration,  or  to  the 
braking  deceleration  which  may  be  considered  as  only  a necessary  experiment- 
al evil,  and  not  a factor  in  the  experiment. 

The  experimental  design  and  resulting  acceleration  levels  at  fracture 
as  well  as  the  spinal  mode  in  which  fracture  occurred  are  listed  in  Table 
1.  A summary  of  average  peak  sled  acceleration  required  to  cause  fracture 
for  all  modes  and  average  age  at  death  of  those  cadavers  is  contained  in 
Table  2.  The  extended  mode  shows  roughly  a 50%  increase  in  the  peak  sled 
acceleration  required  to  cause  fracture  over  the  level  for  the  erect  mode, 
and  a 100%  increase  over  the  flexed  mode.  Yet  the  seat  angle,  the  line  of 
thrust  relative  to  the  pelvis,  and  the  cadaveric  back  structure  itself  were 
all  unchanged,  and  the  instrumentation  and  cadaveric  preparation  were  iden- 
tical in  the  mnioritv  of  the  runs. 

Reference  to  Table  1 reveals  that  the  peak  acceleration  at  fracture  for 
cadaver  11  in  the  extended  mode  is  unknown  since  no  fracture  occurred, 


despite  peak  sled  acceleration  of  24.5  p,  Yet  no  cadaver  run  In  either  the 
erect  or  flexed  mode  escaped  fracture.  More  Importantly,  examination  of 
the  time  of  fracture  (determined  from  the  strain  pages)  shows  that  all  erect 
and  flexed  mode  fractures  occurred  within  the  duration  of  the  sled  experi- 
mental acceleration  pulse,  while  none  of  the  extension  mode  fractures  defin- 
itely occurred  during  this  period,  and  only  one  (cadaver  12)  could  have  oc- 
curred during  experimental  acceleration.  Instead,  examination  of  Table  1 
shows  that  for  cadavers  9 and  10,  run  in  the  extended  mode,  the  vertebral 
fractures  occurred  after  the  end  of  sled  acceleration  and  thus  well  into  the 
period  of  the  braking  deceleration  pulse;  and  therefore,  these  present  a 
different  problem.  The  deceleration  pulse,  which  was  necessarily  used  for 
these  experiments  due  to  location  of  the  accelerator  device,  does  not  occur 
during  operational  ejection  seat  use.  Therefore,  a fracture  due  to  the  de- 
celeration event  is  only  an  experimental  artifact  with  regard  to  the  input 
acceleration  of  interest.  The  deceleration  fractures  that  occurred  to  cad- 
avers 9 and  10  apparently  v/ere  due  to  the  disruption  of  the  anterior  vertebral 
ligament  incidental  to  instrumentation  of  the  anterior  vertebral  bodies  TU- 
LA, which  altered  the  vertebral  column  dvnamic  response  on  deceleration.  The 
time  of  occurrence  of  the  other  fracture  in  the  extended  mode  (cadaver  12)  is 
not  known,  and,  therefore,  has  been  arbitrarily  assigned  to  the  group  due 
to  the  acceleration  pulse  for  conservative  analysis  of  the  data.  The  differ- 
ence between  the  average  fracture  levels  in  the  erect  and  flexed  modes  as 
compared  with  the  extended  mode,  therefore,  could  oossihlv  be  more  signifi- 
cant if  the  reverse  assignment  were  made. 

A t test  was  performed  for  the  fracture  g levels  between  the  various 
spinal  modes  and  the  results  are  given  in  Table  3.  The  differences  in  g 
level  between  the  extended  mode  and  the  other  two  modes  were  found  to  he 
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statistically  significant  (P  = 0.n5).  The  null  hypothesis  vas  rejected 
despite  conservative  estimates  mace  for  the  fracture  level  in  the  extended 
mode. 

When  the  fracture  g levels  in  the  various  modes  were  analyzed,  it 
was  not  possible  to  use  paired  sets  of  data;  so,  the  appropriate  values  of 
t were  obtained  from  the  equation  for  unpaired  data  with  unequal  samples. 

In  particular,  for  extended  mode  the  four  fracture  levels  were  obtained 
under  slightly  different  conditions,  in  that  in  two  of  the  cases  the 
anterior  ligament  was  left  intact  while  it  was  disrupted  in  the  other  two. 

Similarly,  when  the  cadaver  was  not  eviscerated,  the  fracture  level  would 
be  expected  to  decrease.  Therefore,  the  actual  difference  and  the  prob- 
ability that  it  did  not  occur  by  chance  can  only  be  higher  than  those  given 
here. 

Comparison  of  the  fracture  levels  obtained  in  this  series  (Table  2) 

| must  be  made  with  those  reported  by  Ruff  (7).  Table  VI-8  from  his  chapter 

shows  maximum  tolerance  of  the  individual  vertebrae  for  T12-L1  tested 
approximately  in  the  erect  mode  to  be  24.5  g and  23.0  g,  respectively,  and 
minimum  tolerance  for  the  same  vertebrae  as  18.6  g and  18.2  g respectively. 

Yet  the  findings  in  the  present  study  indicate  an  average  tolerance  in  the 

erect  mode  of  10.4  g + 3.79  g in  the  extended  mode  of  17.5  g + 5.55  g.  It  i 

I 

is  believed  that  this  seeming  inconsistency  can  be  resolved  on  the  basis  of 
age.  Ruff's  specimens  were  obtained  at  least  in  part  from  accident  victims. 

i 

It  is  presumed  that  those  accident  victims  were  youthful,  whereas  the  I 

average  age  of  cadavers  in  the  erect  and  extended  modes  of  the  present  ] 

study  was  61  years.  Figure  7 of  the  study  by  McElhaney  and  Roberts  (8), 
shows  that  strength  of  the  vertebral  body  In  the  sixth  decade  of  life  is 
approximately  half  that  in  the  second  decade.  Empirical  data  from  aircraft 
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accidents  indicate  that  the  majority  of  individuals  suffering  vertebral  frac- 
ture are  in  their  twenties. 

If  this  relationship  holds  true,  therefore,  the  average  fracture  level 
for  the  erect  mode  at  age  20,  which  is  the  one  of  interest  extranolnted  from 
the  present  study,  would  be  roughly  20  g-25  g and  for  the  extended  mode 
would  be  35  g-44  g.  Since  cadaveric  bone  is  not  as  strong  as  living  human 
bone  or  fresh  cadaveric  bone,  the  comparison  becomes  potentially  even  more 
meaningful. 

Roentgenograms  of  the  fractures  occurring  to  various  cadaveric  subjects 
showed  anterior  compression  fractures  for  two  cadavers  tested  in  the  flexed 
mode;  for  two  cadavers  in  the  erect  mode.  The  fractures  that  occurred  in 
the  extended  mode,  two  cadavers  were  different  from  those  in  the  erect  or 
flexed  mode.  They  were  somewhat  unexpected  from  previous  pathological  data 
and  were  apparently  due  to  tension  on  the  superior  end  plate  of  the  verte 
bral  body  by  the  attachment  of  the  intervertebral  disk  and/or  ligament  to 
the  end  plate.  If  the  hypothesis  being  tested  were  true,  no  fractures  would 
occur  in  the  posterior  vertebral  structures.  This  was  verified  since  there 
was  no  roentgenographic  evidence  of  damage  to  the  posterior  structures  of 
the  vertebrae  for  any  experimental  acceleration  exposure. 

IV.  THE  ROLE  OF  ARTICULAR  FACETS 

1.  Seat  Pan  Load  for  Different  Spinal  Modes 

Tlie  hyperextension  experiments  not  onlv  verify  the  Ewing  hypothesis, 
hut  also  yielded  important  strain  data  as  shorn  in  Table  4.  In  manv  cases, 
the  decrease  in  anterior  strain  is  a result  of  hyperextending.  The  spine 
is  statistically  significant  at  the  95%  confidence  level.  This  decrease 
is  attributable  to  two  possible  causes.  First,  the  hyperextension  device 
could  interact  with  the  torso  and  ca>-rv  some  of  the  inertial  load  and  thus 
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produce  a drop  in  the  anterior  strain.  Alternately,  the  posterior  structures 
of  the  spine,  principally  the  articular  facets,  could  he  load-hearing  elements 
and  orovide  a second  load  path  for  the  vertebral  column.  The  first  possibil- 
ity was  eliminated  by  comparing  the  seat  pan  loads  of  hvperextended  and 
erect  runs  made  at  the  same  g level.  The  traces  were  identical  as  long  as 
the  run  conditions  were  repeated  exactly.  Thus,  the  hyperextension  block  did 
not  carry  part  of  the  inertial  load,  and  that  the  spinal  load  was 
merely  redistributed.  The  Ewing  hypothesis  could  then  be  extended  to  state 
that  the  articular  facets  constituted  a second  load  path.  The  quantification 
of  this  facet  load  is  described  in  this  section. 

2.  Anatomy  of  the  Facet  Joint 

In  a typical  vertebra,  the  articular  facets  or  processes  are  located  near 
the  junction  of  the  pedicles  with  the  lamina.  The  pair  of  superior  facets 
spring  upward  from  the  pedicles  and  face  in  the  general  posterior  direction, 
while  the  two  inferior  facets  project  downward  from  the  lamina  and  face  anter- 
iorly. The  articular  surfaces  are  lined  with  hyaline  cartilege  and  form  a 
plane  synovial  joint. 

It  is  obvious  from  this  anatomical  arrangement  that  the  overlapping  facets 
perform  the  important  function  of  limiting  rotation  and  of  preventing  one  verte- 
bra from  sliding  with  respect  to  its  adjacent  vertebrae.  The  question  of  whether 
they  are  capable  of  transmitting  compressive  loads  in  the  longitudinal  direc- 
tion of  the  vertebral  column  has  never  really  been  answered.  Ir  most  texts  of 
anatomy,  the  vertebral  body  is  considered  to  be  the  weight-bearing  structure 
of  the  column  (9,  10).  These  references  exemplify  opinions  expressed  in  J 948 
and  1970,  respectively.  The  facets  have  been  said  to  carry  no  load  at  all  (11, 
12).  On  the  other  hand,  Strasser  (13)  and  Nachemson  (14)  have  indicated  that 
the  facets  can  support  a portion  of  the  load  borne  by  the  spine.  However,  in  a 
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later  paper,  Nachemson  (15)  retracted  Ills  earlier  statement  and  declared 
that  the  facets  carry  no  load.  Ills  studies  were  based  on  the  measurement 
of  lntradtscal  pressures  in  isolated  spinal  segments  which  were  subjected 
to  axial  loads  while  the  disk  was  tilted  up  to  5 legrees. 

3.  Qualitative  Evidence  of  Facet  Loads 

Several  techniques  were  used  to  deduce  the  lond-hearinp,  capability  of 
the  articular  facets,  since  direct  measurement  of  force  in  a limited  space 
environment  is  still  beyond  the  state  of  the  art. 

Strain  gages  were  an  effective  means  of  providing  a qualitative  indic- 
ation of  facet  load.  They  were  mounted  on  the  pedicles  and  lamina.  Methods 
of  installing  them  on  vertebral  surfaces  were  developed  to  allow  the  measure- 
ment of  strain  just  about  anywhere  on  the  external  surface  of  a vertebra. 

A detailed  description  of  the  techniques  used  to  install  strain  gages  on 
vertebral  bodies  has  been  described  above.  The  same  methods  were  employed 
for  installing  them  on  the  pedicles,  for  which  an  anterior  approach  is  used. 
To  mount  them  on  the  posterior  surface  of  a lamina,  a posterior  approach 
was  taken,  but  the  hasic  techniques  remained  unchanged. 

As  an  example  of  qualitative  evidence  of  facet  loads,  the  strain  records 
of  two  8 g runs  are  compared.  Figure  9 shoves  the  anterior  (A)  and  posterior 
(P)  strain  gage  output  with  the  spine  in  the  erect  mode.  The  posterior 
gages  were  installed  on  the  lamina  of  TIO,  Til,  and  1,2.  Those  on  T10  and  Til 
were  in  tension  throughout  the  pulse  while  PL2  .tally  in  compression 

and  went  into  tension  at  about  50  ms  after  * iset  of  acceleration.  fata 
from  an  Identical  run  hut  wi th  the  spine  hynorextended  ire  shown  in  figure 
10.  PL2  remained  in  compression  for  the  entire  duration  of  the  pulse, 
while  PT10  and  PT11  underwent  a change  in  sign  from  compression  to  tension. 

A marked  reduction  in  anterior  strain  is  also  nuite  evident. 
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Fig.  10  Strain  data  for  a run  in  ;nc^ 
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4.  Ouanti tative  Evidence  of  facet  Loads 

Ouantitative  measures  of  facet  load  were  obtained  by  means  of  an  in- 
tervertebral load  cell  (IVLC)  designed  to  fit  under  a disk  of  a vertebral 
body.  it  can  measure  both  axial  force  and  the  eccentricity  of  that  force 
with  respect  to  its  geometric  center.  The  first  model  of  an  IVLC  was  almost 

I 

1 in.  thick  and  is  believed  to  be  the  first  transducer  capable  of  actually 
measuring  the  load  carried  by  a vertebral  bodv  in  an  intact  vertebral 
column.  A thinner  version  is  shown  in  Figure  11.  This  load  cell  was  de- 
signed to  fit  above  a lumbar  disk,  replacing  the  inferior  segment  of  a 
vertebral  body.  A douhle-bladed  rotary  saw  was  constructed  to  cut  slots 
of  a precise  width  across  the  spine  to  accomodate  the  IVLC  without  damaging 
the  neural  arch.  The  lin  shown  in  Figure  11,  enables  a strap  to  hold  the 
IVLC  in  place  during  the  experimental  run.  A properly  installed  IVLC  will 
not  result  in  a change  in  length  of  the  ve.  -ebral  column,  and,  if  the  thin- 
ner model  is  used,  the  mobility  of  the  column  is  also  not  affected.  When 
an  IVLC  is  used,  strain  gages  are  installed  on  both  the  anterior  and  posterior. 

The  purpose  of  the  IVLC  was  to  measure  the  load  carried  by  the  verte- 
bral body  and  to  compare  it  with  that  borne  by  the  column,  the  total  spine 
load.  The  latter  was  taken  to  be  proportional  to  the  measured  seat  pan 
load.  The  ratio  used  was  the  weight  of  the  torso  above  the  IVLC  to  the 
total  body  weight.  Any  difference  between  the  spine  load  and  the  IVLC 
output  is  the  load  carried  by  the  articular  facets,  the  facet  load.  Justi- 
fications for  the  validity  of  this  method  of  deducing  facet  load  are  pro- 
vided in  the  discussion  section  of  the  paper, 

IVLC  data  were  obtained  from  three  different  cadavers.  The  initial 
thicker  version  was  used  In  the  first  two  cadavers,  while  tiie  second  thinner 
model  was  employed  during  runs  made  on  the  third. 
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The  various  spinal  loads  for  a 10  g run  in  the  erect  mode  are  shown  in 
Figure  12.  The  IVLC  replaced  the  inferior  portion  of  the  L3-L4  disk  and 
the  sunerior  segment  of  L4.  The  neural  arch  of  T.4  was  left  intact.  The 
facets  were  in  compression  for  the  first  125  ms  of  the  pulse.  However, 
as  the  head  and  torso  rotated  forward,  the  facets  unloaded  and  went  into 
tension,  resulting  in  an  intervertebral  body  force  larger  than  the  total 
spine  load.  The  facet  load  and  PL3  strain  also  show  good  correlation, 
with  the  zero  crossover  point  of  both  traces  occurring  almost  simultaneously. 

A series  of  runs  were  carried  out  in  a cadaver  using  the  thinner  model 
of  the  IVLC.  It  was  placed  about  the  L3-L4  disk  by  replacing  the  inferior 
segment  of  L3.  The  disk  was  virtually  intact.  Figure  13  shows  data  from 
two  8 g runs  made  in  the  erect  and  hyperextended  mode.  In  the  first  case 
the  facets  unloaded  rather  early  in  the  acceleration  pulse  (at  about  60  ms) 
but  there  was  confirmation  from  PL3  strain.  In-  the  hyperextended  mode, 
the  facets  remained  in  compression  for  the  entire  pulse  duration,  as  shown 
in  Figure  13.  PL3  and  PL4  strains  were  also  in  compression  throughout  the 
pulse.  The  total  spine  load  was  the  same  in  these  two  runs,  but  the  inter- 
vertebral body  force  decreased  by  about  400  lb  as  a result  of  hyperextension. 
These  runs  were  repeated  at  6 and  10  g and  the  results  were  essentially 
the  same  as  those  at  8 g. 

In  summary,  the  following  conclusions  can  be  listed: 

1.  The  articular  facets  are  capable  of  bearing  compressive  and  tensile 
loads.  When  the  facets  are  in  tension,  it  is  assumed  that  the  joints  and 
the  associated  capsular  attachments  are  being  stretched. 

2.  Strain  gages  were  employed  to  provide  qualitative  evidence  of  facet 

load. 

3.  Intervertebral  body  force  was  measured  in  an  intact  spine  during 
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Fig.  12  Intervertebral  and  facet 
loads  for  an  erect  mode  run 


TIME  (msec) 


£ 1 1 
2 1 
£ ip 

u iuj 

oc  , cc 

LJ  ' LU 

t I ~ 

x»  >! 

;i  \ 


UJ 

H ^ 
X h* 

£ G 

c 111 

111  OC 

£ “ 

1 ^ 
, IO 

S s! 

2 n 


UJ 

h 

X 

UJ 

-nr 

Q-  UJ 

^ CL 
>- 


CO 

TJ 

CO 

O <b 

*-)  CO 
■U  <D 

0)  CD 
U *H  O 
CO  B 

4-|  OJ 
TO  TO 
4-1  CO  CD 

o Btj 
c 

C l/l  II 
O C 
CO  3 x 
*H  M QJ 
Li  Li 
CO  Li  CD 
CL  O CL 
S 4-i  >, 
O JS 

o c 

CO  G 
on  Li  CO 


r , , 

f , , , 

1 » 1 

, . r , 

fa  « 

1 1 

0 « 0 

OOOO 
OOO 
cvj  it  co 

1 1 1 

0 

-250 

-500 

-750 

0000 

8S  S 

1 

O O O 

0 0 

lO  id 

, 1 

( *)  Nivdis 

(6)  "1300V 

(Ql)  30U0J 

(qi)  aouod 

(QD30U03 

I 


\ 

impact  acceleration  by  means  of  a specia]lv  designed  intervertebral  load 
cell  (IVLC) . 

4.  From  the  IVLC  and  seat  pan  load  cell  output,  a facet  load  historv 
was  computed. 

5.  A better  understanding  of  the  injury  mechanism  of  the  snine  has 
been  achieved. 

I 

6.  Hyperextension  of  the  spine  transfers  more  load  to  the  facets. 

7.  The  proportion  of  the  load  carried  bv  the  facets  appear  to  increase 
with  the  decreasing  g levels,  suggesting  that  they  mav  also  carry  a portion 
of  the  static  body  weight  when  the  bcdv  assumes  normal  erect  posture. 
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Table  1 Fracture  Levels  and  Spinal  Modes 


Cadaver 

No. 

Rate  of 
Onset 
(g/sec) 

Peak 

AccL 

(8) 

Duration 

(m/sec) 

Time  After  w . 
n . Mode  At 

Onset  „ 

, , . Fracture 

(m/sec) 

No 

Modes 

. of 
Tested 

Fractured  Vertebra 

1 

246 

11.0 

195 

116 

erect 

2 

T12 

2 

320 

7.5 

238 

138 

erect 

2 

Tl  0 

3 

375 

5.5 

300 

130 

erect 

2 

TH 

4 

Data  Lost 

14.0 

Data  Lost 

110 

erect 

2 

T9 

5 

418 

14.0 

208 

60-140 

erect 

2 

Til 

6 

830 

7.0 

263 

120 

flexed 

3 

L2 

7 

162 

11.0 

248 

85 

flexed 

3 

Til 

8 

333 

9.0 

245 

195 

flexed 

3 

L2M.4 

9 

267 

14.0 

206 

390 

extended 

1 Tension  fracture  on 

sunerior  end  nlate 
of  LI 

10 

483 

20.0 

170 

290 

extended 

1 Tension  fracture  on 

sunerior  end  plate 
of  Tl 2 5,  L2 

11 

280 

24.5 

140  no  fracture  extended 

1 

No  fracture 

12 

333 

12.5 

178 

not  known 

extended 

1 

T8 

26 
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Table  2 Summary  of  Peak  Acceleration  Values  at  Fracture  in  the  Three 

Spinal  Modes 


Fracture 

Level 

(g) 


No . of 
Cadavers 


Average 

Age 

(vrs) 


Extended 

Erect 

Flexed 


17.75  + 5.55 
10.4  + 3.79 
9.0  + 2.00 


Table  3 Student's  t Test  of  Fracture  g-Levels  Between  the  Spinal  Modes 


Modes 


Sample 

Size 


Table  4 Paired  t Test  of  Strain  Data 


Level  n 


Vertebra 
AL1  AL2 

t P n t 


4.5  5 1.63  0.18  7 2.29  0.06  5 3.63  0.02  7 8.03  >.01  5 2.88  0.04  5 3.87  0.02 

7.5  3 4.48  0.05  3 6.49  0.02 

9.0  3 5.86  0.03  3 19.22  >.01  3 6.97  0.04  3 5.74  0.03  3 3.08  0.09  
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V.  THE  EFFECT  OF  ABDOMINAL  PRESSURE 

The  ability  of  the  abdominal  cavity  to  act  as  a third  load  nath  was 
investigated.  Its  role  as  a load-bearing  pressurized  cylinder  during 
weight-lifting  has  been  documented  by  Morris  et  al  (16). 

Three  male  cadavers  were  subjected  to  +0  acceleration  in  a seated 
position.  A detailed  description  of  the  methodologv  has  been  gi''en  by 
Prasad  et  al  (17).  An  intervertebral  load  cell  (IVLC)  was  positioned  in  a 
slotted  lumbar  vertebral  body  at  either  1,2  or  L3.  It  measures  the  load 
carried  by  the  body  and  its  line  of  action.  A rubber  balloon  containing 
a miniature  pressure  transducer  was  placed  in  the  eviscerated  abdominal 
cavity.  The  input  pulse  was  trapezoidal  ir  shape  and  the  peak  acceleration 
ranged  from  A to  9 g.  Its  duration  was  about  200  ms.  The  effect  of  abdomin- 
al pressure  on  spinal  load  distribution  was  studied  by  making  runs  with  and 
without  pre-pressurization  of  the  balloon.  The  initial  static  pressures 
varied  from  6.9  to  20.7  kPa.  During  the  runs,  seat  acceleration,  seat-pan 
load,  intervertebral  load  and  abdominal  pressure  were  measured. 

To  determine  the  facet  force^^,  the  vertical  forces  acting  on  the 
torso  at  the  IVLC  level  were  equated  to  the  inertial  force  (IF) : 

RB  + RF  + FB  - RM  - W = IF  (1) 

where  RB  is  the  measured  intervertebral  force  and  the  inertial  force  was 
taken  to  be  proportional  to  the  seat-pan  load.  It  was  reduced  by  the  ratio 
of  the  body  weight  above  the  IVLC  (W) , to  the  total  body  weight.  The 
abdominal  force  (FB) , was  computed  as  the  product  of  the  abdominal  pressure 
and  the  abdominal  cross-sectional  area  at  the  level  of  the  IVLC.  The 
abdominal  wall  reaction  (RM) , was  assumed  to  be  pronortional  to  the  abdominal 
force.  The  proportionality  factor  was  allowed  to  vary  from  0.15  to  0.50 
without  affecting  the  facet  force  apprec'ably.  A value  of  0.3  gave  reasonable 
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estimates  of  RM  which  were  less  than  200  N.  The  facet  force  computed  from 
Equation  1 is  plotted  in  Figure  1 h along  with  the  intervertebral  load.  They 
are  for  a 9 g run  with  and  without  nre-pressurization  of  13.8  pKa.  The 
facets  initially  shared  in  supporting  the  inertial  load.  However,  as  the 
head  and  torso  rotated  forward,  they  unloaded  and  went  into  tension.  The 
effect  of  the  abdominal  load  path  was  to  reduce  the  amount  of  facet  unload- 
ing by  developing  a moment  against  forward  rotation.  The  intervertebral 
force  was  also  reduced  by  this  third  load  path.  The  results  were  similar 
in  the  runs  made  on  the  other  two  cadavers. 

Although  abdominal  pressure  altered  the  load  distribution  across  the 

spine,  the  basic  dynamic  response  described  by  Prasad  et  al  (17)  was  not 

changed.  An  abdominal  pressure  of  13.8  pKa  (104  mm  of  Hg)  is  difficult  to 

generate  in  200  ms  but  an  acceleration  level  of  9 g is  less  than  half  of 

that  necessary  for  pilot  ejection.  Thus;  the  effect  of  abdominal  pressure 

is  only  of  significance  at  low  g leve.ls,  and  cannot  affect  the  loads 

sustained  by  the  spine  during  ejection.  The  computation  of  facet  force 

did  not  account  for  the  balance  of  angular  momentum  due  to  the  lack  of 

a. gular  acceleration  data  for  the  head  and  torso.  An  estimate  of  EM  was 

made  to  eliminate  the  second  unknown  in  Equation  1. 

VI.  MUSCHLAR  RESPONSE  TO  +G  ACCELERATION 

z 

1.  Response  of  the  Living  System 

The  biodynamic  response  of  the  living  spine  to  +C>z  impact  acceleration 
has  been  studied  experimentally  by  many  investigators.  They  have  used  a 
variety  of  test  subjects,  including  human  volunteers,  cadavers,  and  animals, 
ranging  from  brown  bears  to  rhesus  monkeys.  Early  work  bv  Geerts  (18),  Love- 
lace et  al  (19),  and  Stapp  (20,  21)  was  carried  out  with  human  volunteers  to 

establish  human  tolerance  limits- to  +G  acceleration. 

z 


4 


___ . 





Cadaver  2627 
9g.Run,  IVLC  in  L2 


Fig.  14  Spinal  loads  with  and 
without  abdominal  pressure 


30 


Animal  studies  were  carried  out  by  Beeding  and  Cook  (22),  who  used 
swine  and  bears,  and  by  Stapp  (23),  who  tested  swine  and  chimpanzees. 

The  principal  aim  was  to  establish  injury  patterns  above  the  tolerance 
limits  and  the  animals  were  all  anesthetized.  Recently,  Kazarian  et  al  (24, 
25) , reported  the  biodynamic  response  of  the  rhesus  monkev  and  discussed 
spinal  fracture  patterns  sustained  during  high  amplitude  short  duration 
+Cz  impact.  The  study  was  extended  by  Kazarian  (26)  to  include  data  relat- 
ing to  chimpanzees  and  baboons. 

In  this  study,  cadaveric  responses  of  the  spine  have  been  carefully 
documented.  To  complete  the  experimental  phase,  it  was  necessarv  to  in- 
vestigate spinal  muscular  response  to  -l-C^  impact  acceleration.  In  part- 
icular, the  time  delay  of  muscle  activation  and  the  form  of  muscle  force- 
time curve  are  reported.  Electromyographic  signals  from  canine  musculature 
were  acquired  during  low  level  +0^  acceleration.*  Spinal  electromyographic 
activity  was  anticipated  as  a result  of  stretching  of  spinal  muscles,  thus 
bringing  about  muscular  contraction. 

2.  Development  of  Instrumentation 

The  experiments  were  carried  out  in  the  same  vertical  accelerator. 

The  sled  was  modified  to  accommodate  a dog,  which  was  held  in  a "seated" 
position  by  a combination  of  belt  restraints  and  padded  side  boards.  A 
specially  designed  seat  back  contained  a window  through  which  electro- 
myographic leads  could  be  connected  to  the  preamplifiers.  Trailing  cable 
conveyed  the  signals  to  a data  recording  station  located  on  the  fourth 
floor. 

* The  experiments  reported  in  this  paper  were  conducted  in  accordance  with 
the  "Guide  for  the  Care  and  tlse  of  Laboratory  Animals,"  DREW  Publications  No. 
(NIH)  74-23,  Revised  1972. 
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One  of  the  major  difficulties  in  recording  noise-free  electromyographic 
signals  from  an  impact  sled  is  the  availability  of  shock  resistant  on-bonnl 
amplifiers.  They  are  required  to  increase  the  s I gua  I-  t o tin  i ■•,<•  ratio  [it  Ini 
to  the  transmission  oi  the.  electromyographic  signals  througii  20  ro  of  trail 
ing  cables.  Such  amplifiers  were  difficult  to  obtain  commercially  and  a 
set  of  eight  was  designed  and  fabricated  in  the  laboratory.  As  shown  in 
Figure  15,  tliev  were  packaged  within  a single  case,  which  could  l.e  mounted 
on  the  sled.  These  were  differential  amplifiers  with  a variable  gain  from 

0 to  about  10,000,  and  the  output  was  filtered  by  second  order  active  band 
pass  filters  with  both  low  pass  and  high  pass  provisions.  The  amplifier 
cut  off  frequencies  above  600  Hz,  and  below  200  llz,  since  the  frequency 
content  of  canine  e lectromyographic  signals  range  from  200  to  400  Hz.  The 
filters  and  the  differential  amplification  used  eliminated  most  of  the  noise. 
The  latter  feature  also  permitted  freedom  of  choice  in  the  placing  of  the 
electrodes  to  record  only  the  difference  in  potential  between  them.  The 
tendency  for  the  amplifier  to  become  saturated  was  curbed  by  placing  a 

1 M fl  resistor  across  the  input.  The  amplifiers  were  checked  out  by 
acquiring  electromyographic  data  from  dogs  walking  on  a treadmill.  The 
wire  electrodes  were  inserted  in  leg  muscle  groups  known  to  be  active 
during,  walking.  The  electrode  material  and  insertion  techniques  were 
similar  to  those  developed  by  Basmajian  (27).  By  visually  correlating 
the  oscillograph  records  with  the  gait  of  the  dog,  it  was  possible  to 
distinguish  the  electromvogrnphi c signals  from  noise.  From  these  tests 
the  necessary  adjustments  and  modifications  were  made  to  maximize  the 
signal-to-noise  ratio. 

f'eedle  electrodes  inserted  at  a shallow  angle  into  (be  muscle  were 
found  to  tie  a stable  in;t  illation.  The  needle  shank  was  boni  10  to  15  deg, 
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rees  prior  to  Insertion  to  permit  deeper  penetration,  and  a circumferen- 
tial groove  was  machined  Into  the  handle  of  the  needle  so  that  it  could 
be  securely  tied  in  place  hv a suture. 

The  unipolar  electrode  was  selected  for  the  dynamic  exper iments . * 

It  was  37  mm  long  and  had  a tapered  stainless  step1  cannula,  wit*  a co- 
axial insulated  nickel-chrome  alloy  wire  exposed  at  the  tip.  Mien  com- 
pared wi tli  a similar  bipolar  electrode,  the  signals  were  of  merely  identi- 
cal wave  form  when  isometric  contractions  of  a human  flexor  carpi  radial  is 
muscle  were  recorded.  However,  dual  monopolar  electrodes  gave  signals 
that  were  four  times  as  strong  as  those  from  a bipolar  electrode. 

An  excellent  shielding  and  grounding  system  was  required  to  minimize 
noise  in  the  minute  electromyographic  signals.  Each  electrode  was  shield- 
ed by  the  needle  cannula  and  a shielded  coaxial  cable  was  used  for  each 
needle.  Ground  loops  were  carefully  avoided.  For  each  channel  the  two 
active  lead  cable  shields,  the  body  reference  ground,  and  the  direct  cur- 
rent power  supply  ground  were  connected  in  parallel  to  the  electromyo- 
graphic preamplifier  reference  ground.  This  ground  was  carried  t>  rough 
the  trailing  cable  shield  to  the  common  instrumentation  ground  at  the 
data  recording  station.  All  ground  and  shields  were  open  ended  ind  were 
carried  to  a common  ground  sink. 

This  impact  acceleration  experiment,  had  an  additional  noise  problem, 
which  was  traced  to  the  whipping  of  the  electrode  leads,  especially  If  they 
struck  a rigid  surface.  Moise  levels  roughlv  equivalent  to  those  of  the 
electromyogram  were  produced.  To  minimize  this  noise,  several  tvpes  of 
miniature  coaxial  cables  were  tested.  It  was  determined  that  type  RG-l'*6A/u 
produced  the  least  noise  during  the  run.  Tt  was  also  essential  that  each 


* Manufactured  by  Teca  Corporation,  White  Plains,  Nevt  Yorl  10A03. 


cable  be  securely  tied  down  to  all  surfaces  that  it  could  hit. 


3.  Electrode  Sites  and  Experimental  Puns 

The  muscle  groups  chosen  for  electromyographic  study  were  known 
to  be  active  during  flexion  and  extension  movements  of  the  spine.  The 
following  groups  were  monitored : 

1.  Erector  Spinae 

a.  Tliocostalis  (lumbar) 

b.  Lonf,issimus  (lumbar,  thoracic,  and  cervical) 

c.  Spinalis  (thoracic,  cervical) 

2.  Multifudus  (lumbar,  thoracic) 

3.  Semi  spinalis  (cervical) 

Approximately  24  check-out  runs  were  made  to  test  various  components  of  the 
experimental  system  and  to  develop  techniques  for  acquiring  dynamic  electromyo- 
graphic signals.  The  on-board  preamplifiers  were  subjected  to  runs  at  5 g to 
verify  their  structural  integrity  and  to  ensure  that  no  spurious  signals  were  gen 
erated  as  a result  of  the  shock  loading.  There  were  also  several  runs  desig- 
nated to  eliminate  ground  loops  and  other  sources  of  noise.  Four  runs  were  then 
made  with  the  entire  electromyography  system  on  board,  but  the  electrodes  were 
not  inserted  into  any  muscle;  thev  were  kept  about  1 cm  apart  in  air.  The  noise 
level  was  found  to  be  low  and  there  were  no  spurious  signals. 

A total  of  36  runs  at  3 and  4 g were  made  with  two  additional  canine 
subjects  to  acquire  electromyographic  data.  They  were  male  beagles  weigh- 
ing approximately  10  kg  each.  The  experimental  procedure  consisted  of 
setting  up  the  instrumentation  to  record  up  to  eight  channels  for  electro- 
myography, the  sled  acceleration,  and  the  seat-pan  axial  load  and  sagittal 
bending  moment.  Lateral  supports  were  provided  in  addition  to  a cross  chest 
belt,  and  upper  torso  Y-yoke  restraint,  and  a lap  belt. 

Upon  completion  of  a series  of  runs,  the  dog  was  given  0.5  ml  of 
Acepromazine  and  postrun  x-ray  films  were  taken  to  verify  needle  electrode 
positions.  The  animal  was  also  given  1 ml  of  Floricilin  intramuscularly 
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after  the  electrodes  and  sutures  were  removed. 

Figure  16  shows  typical  electromyographic  data  from  the  multifudus 
1 umbo rum.  The  input  acceleration  pulse  is  superimposed  to  show  the  del  ,<v 
in  muscle  activation.  The  plot  also  shows  that  prior  to  the  run  there 
was  very  little  electromyographic  activity. 

4.  Data  Analysis  and  the  Intergrated  EMC 

The  analog  data  were  digitized  at  a rate  of  2000  samples  per  second  per 
cnannel,  using  a PDP-8E  minicomputer.  The  data  were  prefiltered  at  a 10oq  Pz 
level  to  avoid  "aliasing"  or  the  generation  of  false  low  frequenev  components. 
Hannan  et  al  (28)  and  Viitasalo  and  Komi  (29)  have  reported  sampling  rates  of 
1000  and  1600  per  second,  respectively.  The  digitized  data  were  recorded  on 
magnetic  tape  by  an  industry-compatible  digital  recorder  so  that  thev  could  be 
processed  by  an  IBM  360/67  computer. 

"intergrated"  electromyography  has  long  been  recognized  as  a measure  of 
muscle  force.  The  proportionality  factor  has  been  difficult  to  establish  except 
for  special  cases  of  isometric  contraction  or  isotonic  contraction  at  a constant 
velocity.  To  obtain  this  force,  the  electromvographic  signal  is  rectified  and  a 
curve  smoothing  procedure  is  applied  to  the  rectified  data.  In  this  study  the 
rectification  was  performed  by  means  of  a root-mean  square  (RMS)  analysis,  based 
on  the  frequency  content  of  the  signal.  For  an  analog  alternating  signal,  its 
RMS  value  is  given  by, 

E=  [(1/T)  /T  e2dt]?i 
o 

where  E=RMS 

e“Instantaneous  value 

T=Time  period  of  the  wave  form  window 

The  equivalent  digital  form  is 

E-  [(1/N)  r e 2 d 1 1 * ’ 

*-  i-1  1 

where  N=number  of  points  in  the  window. 


Hv  using  » window  size  of  5(N*-5)  , the  RMS  analysis  covered  at  least  one 
fourth  of  the  period  of  a 100  Hz  waveform,  • i rice  the  di  gi  t i y.i  rig  rate  was 
2000  per  second.  The  RMS  data  were  then  filtered  hv  the  ITT  method,  with 
•i  low  pass  cut-off  frequency  of  20  ||/.  I'hls  ensured  that  the  natural 
frequency  of  8 to  12  Hz  for  physiological  tremor  in  muscle  was  preserved. 
fHntmal  phase  shift  resulted  from  the  u.se  of  this  digital  filter. 

Figure  17  and  18  are  examples  of  electromvographical.lv  derived  force 
for  each  of  the  nine  muscle  sites  monitored  in  the  two  canine  suhlects. 

In  general,  there  was  a major  burst  of  electromyographic  activity  short! 
after  the  onset  of  acceleration.  This  initial  hurst  was  usually  the 
strongest  and  resulted  in  the  development  of  a peak  electromvographical 1 v 
derived  force.  A sustained  electromyographic  signal  and  thus  a sustained 
force  is  shown  in  Figure  18;  these  were,  rather  uncommon  and  appeared  to  be 
associated  with  prerun  muscle  tone. 

In  terms  of  relative  magnitudes,  the  iliocostalls  lumlorum  had  the 
highest  peak  value,  at  both  3 g and  5 g.  Table  5 lists  mean  peak  values 
for  both  dogs  at  two  g levels.  On  the  average,  the  peak  force  was  higher 
at  5 g than  at  3 g. 

Table  5 Peak  Values  of  Klectromyograpli  i cal  1 v Deri  ved  Forcp 


Muscle  Relative  Peak  Value  d-  S.  D.  (mv) 


r 

r.roup 

3 

R 

Runs 

(N) 

5 

R 

Runs 

(N) 

All 

Runs 

(M) 

Semispin.  cap. 

4 

d- 

2.0 

(8) 

6 

+ 

4.5 

(9) 

5 

d- 

3.6 

(17) 

Spinalis  cerv. 

8 

+ 

9.0 

(14) 

12 

+ 

14.2 

(15) 

10 

+ 

12.0 

(29) 

i 

Spinalis  thor. 

6 

+ 

4 . 4 

(14) 

12 

+ 

10.1 

05) 

o 

+ 

8.2. 

(29) 

* 

bong.  cerv. 

8 

+ 

9.8 

(9) 

7 

+ 

7.0 

(14) 

8 

d 

fi.C 

(23) 

k 

Long.  thor. 

13 

(1) 

1 

0) 

7 

d 

8.6 

(2) 

Long.  lum. 

4 

+ 

3.4 

(20) 

10 

d- 

11.3 

(17) 

7 

+ 

8.4 

(37) 

Mult.  thor. 

5 

d- 

4.8 

(9) 

6 

+ 

4.5 

(V 

6 

4- 

4.6 

(18) 

Mult.  lum. 

4 

+ 

2.1 

(12) 

6 

d 

4.0 

05) 

5 

f 

3.4 

(27) 

1 1 locos t.  lum. 

20 

+ 

19.3 

(2) 

15 

d 

5.1 

(3) 

17 

+ 

10.6 

(5) 
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5.  Muscle  Response  Delay 

Tlie  proprioceptive  response  of  muscles  has  been  studied  by  many  in- 
vestigators. Neural  delay  time  in  the  human  arm  muscle  was  found  to  be 
60  to  80  ms  by  Hammond  (30).  This  delay  includes  the  time  required  to 
stretch  the  muscle  spindles  that  form  the  completion  of  the  reflex  arc. 
Soechting  et  al  (31)  have  also  found  a delay  time  of  90  to  140  ms  in  the 
human  arm,  and  Foust  et  al  (32)  found  reflex  times  for  human  neck  muscles 
to  be  between  54  and  92  ms  from  the  onset  of  head  acceleration.  For  whole 


body  +Gz  acceleration  in  the  dog,  muscular  response  delay  is  defined 
as  the  time  from  onset  of  sled  acceleration  to  the  initial  burst  of  elec- 
tromyographic activity.  Table  6 shows  the  average  delay  time  for  each 
of  the  muscles  monitored.  The  delay  is  shorter  than  those 

Table  6 Average  Electromyographic  Signal  Onset  Delay 


Muscle 
Group 

Semispin.  cap. 
Spinalis  cerv. 
Spinalis  thor. 
Long.  cerv. 
Long.  thor. 
Long.  lum. 
Mult.  thor. 
Mult.  lum. 
Iliocost.  lum. 


Delay  Time  + S.  D.  (ms) 


3 

Runs 

(N) 

5 

P, 

Runs 

(N) 

All 

Runs 

(N) 

25 

+ 

5.7 

(7) 

28 

+ 

12.9 

(9) 

26 

+ 

10.2 

(16) 

21 

+ 

6.8 

(12) 

27 

+ 

11.3 

(14) 

24 

+ 

9.9 

(26) 

32 

+ 

8.8 

(13) 

28 

+ 

9.2 

(14) 

30 

+ 

9.1 

(27) 

26 

+ 

4.2 

(7) 

36 

+ 

10.6 

(11) 

32 

+ 

9.7 

(18) 

— 

(0) 

22 

(1) 

22 

(1) 

36 

+ 

22.3 

(18) 

25 

+ 

8.8 

(17) 

31 

+ 

17.8 

(35) 

37 

+ 

11.7 

(ID 

27 

+ 

12.1 

(8) 

32 

+ 

12.5 

(19) 

58 

+ 

29.4 

(6) 

37 

+ 

12.3 

(ID 

44 

+ 

21.7 

(17) 

28 

(1) 

44 

+ 

41.5 

(3) 

40 

+ 

34.9 

(4) 

of  human  muscles.  The  difference  is  also  attributable  to  the  close  prox- 
imity of  the  spinal  musculature  to  the  cord,  so  that  the  reflex  is  shorter. 

The  large  variation  in  delay  times  in  a given  muscle  grouo  can  be 
due  to  several  factors.  Basmajian  (27)  found  that  conduction  velocities  varied 
widely  from  time  to  time  in  any  individual.  Furthermore , y efferents  from 
the  central  nervous  system  provide  a variable  mechanical  bias  to  the 
spindles  and  can  alter  their  sensitivity  to  stretch.  Finally,  the  spinal 
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configuration  of  the  subject  can  Iwvn  a s * gni  f i v ant  o f rec * on  !’:<■  re:'  n . 
time.  Increased  kvphos  1 s can  shorten  the  dolav  in  the  nul  ti  fudii  ; thoraci' 
because  the  curved  spine  is  compressed  during  +(1  acceleration. 

f) . Conclusions 

The  following  conclusions  were  m .dc  • 

1.  Noise-free  electromvogranh  i c ; ign  : 1 s ''ere  acquired  lu*"ir.  • wlx.tr 
liody  -H^imnact  acceleration.  using  canine  subjects.  Nino  muscle  froims 
of  the  spinal  musculature,  were  monitored. 

2.  The  validitv  ol  the  elect romvographfc  data  was  demonstrated  1 '• 
eliminating  possible  causes  of  noise  and  spurious  signals.  Small  doses 
of  a tranquilizer  administered  two  to  two  and  one-ha1 f hours  prior  to 
the  experimental  runs  did  not  affect  muscular  response. 

1.  The  fast  Fourier  transform  method  of  spectral  ana  I "si-  was  ennlced 
to  identify  electromyographic  signals  and  to  filter  the  root-mean  square 
"intergrated"  electromvogram  to  provide  an  electron"’ogranhlcall v derived 
force-time  curve. 

4.  The  relative  peak  values  of  this  force  were  compared  for  runs  at 
two  different  g levels.  The  delay  times  in  muscular  response  and  for  tire 
electromvographlcallv  derived  force  to  roach  a peal-  value  err o memt  f f 1 < 1 . 
The  delay  in  the  onset  of  Increased  electromvogranhic  activity  ’'uried 
from  21  to  58  ms.  These  values  are  somewhar  less  than  those  observed  In 
human  cervical  muscles.  ,,’1'e  time  for  the  elect  rom,rogranI  Ini'"  •or  ve>* 
force  to  reach  its  pea!  ranged  from  45  to  155  ms. 

VII.  A Th'O-DlMI'NS IONAI.  "A'I'i"" 'ATICAI.  MOnri. 

1.  Introduction 

After  Prasad  et  al  (17)  reporter!  the  existence  of  two  load  paths 
along  the  spine,  one  through  the  vertebral  bodv  and  the  other  throng1'  t hi 
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r.a,  v <i  the  art  i ii'iUtr  facets,  previous  models  are  inadequate  in  the 
co-rlete  rsf ••'illation  of  spin  I rrsnonse.  Hence , a two-dimensional  discrete 
"arn  cter  model  of  tin.  spine  <<as  formula  eel  to  incorporate  the  two  load 
1 at!:-.-  -a  factor  not  considered  in  any  of  the  existing  mathematical  models, 
’.'hi  "Meciive  is  to  propose  a realistic  spinal  model  simulating  dynamic 
rcvvmse  ••Mel  can  he  rot  pared  to  that  of  an  intact  spine  tested  in  vivo. 
Such  a model  should  consider  the  natural  spinal  curvatures  and  the  effects 
of  I'Hxior  and  ccentric  irertial  loading  on  i:  e spine.  Head  and  neck  mo- 
tions are  simulated  and  their  effects  on  the  forces  and  moments  in  the 
thoracic  and  lumbar  spine  can  he  studies  for  off-axis  impacts  in  the  mid- 
sm  : ttal  plane.  The  input  acceleration  pulse  can  he  an  arbitrarv  function 
of  time  and  the  restraint  and  support  systems  have  been  included  to  proper- 
ly simulate  a seated  vehicle  occupant.  In  view  of  the  fact  that  the  experi- 
mental data  for  validation  of  the  model  were  obtained  from  cadaveric  runs 
with  the  spine  in  the  erect  and  hvpei'extended  mode,  the  model  incorporates 
the  ability  to  simulate  both  sninal  configurations. 

2.  Assumptions 

"’he  following  assumptions  were  made  in  the  mathematical  development: 

II.  the  114  vertebral  bodies,  the  head  and  the  pelvis  are  rigid  bodies 
constrained  to  trove  in  the  mid-sagittal  plane. 

2.  Each  rigid  hodv  lias  three  degrees-of-freedom  in  the  mid-sagittal 
plane — two  translational  and  one  rotational . 

'3.  The  intervertebral  disks  ire  massless  and  deformation  of  the 
spine  takes  nlace  at  the  disks. 

4.  The  .lists  are  r. •planed  by  a svstem  of  springs  ami  dampers — one 
spring  and  damp. a-  for  axial  forces,  one  spring,  and  damper  for  shear  forces, 
and  another  spring  and  damper  arrangement  for  restoring  torques  due  to  re- 
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lative  angular  motion  between  adjacent  vrrtehr.il  hosier.  J 

. 

“>.  The  facets  and  laminae  are  springs  connected  t<  the  vertebral  . 

j 

body  hv  a massless  rigid  rod. 

4.  inch  r i >t  i d l-o.lv  i • assumed  to  rnrr"  i oori  ion  - 1 1 1 • . t o 

I 

weight  which  is  eccentric  with  respect  to  the  centerline  of  the  spine. 

7.  The  rigid  bodies  're  arranged  to  simulate  the  « i nil  curvet  nr 
as  closely  as  possible. 

1.  Equations  of  Motion 

Figure  19  shows  the  configuration  of  two  successive  rigid  l:nks  t'  >t 

■i 

have  undergone  relative  translation  and  rotation  due-  to  deformation  o'  the 
disk.  Tn  the  undeformed  state,  it  has  been  assumed  that  the  axis  o'  a.iv 
disk  is  coincident  wi  tit  the  axis  of  the  vertebra,  t'at  is,  the  rigid  Hid 
immediately  below  it.  The  position  in  the  mid-sagittal  plane  of  t die  ent<  r 
of  the  ith  rigid  link  is  determined  1 >y  tlie  three  generalized  coordinates, 
iw  , w^  , 9j  as  sliown  in  the  figure. 

The  reaction  forces  and  moments  developed  on  the  ith  rigid  ’ odv 
due  to  deformation  of  the  ith  and  the  i + I disks  and  the  laminae  are 
shown  in  Figure  ?<).  The  changes  in  the  relative  angular  orientations  u 1 
the  lengths  AiJ,  i*C,  AN!)  AjA^  sliown  in  Figure  19  are  used  to  compute  the 
restoring  moments,  axial  and  shear  forces  on  the  vertebral  1-odv  an!  reac- 
tion forces  at  the  articular  facets.  These  moments  and  forces  are.  ri"nTi 
by  the  following  equations: 


T7Y 

- YKj (AC{ 

“ Ac0i> 

+ CiACi 

(2) 

T7V  j 

= xk  . (nr , 

si  i 

" ^Oi> 

+ CslR‘Ci 

(3) 

Bi  * 

XKT1{(01  - 

V 

(ni-i  - W]  + 'ti((:i 

*l-i> 

(4) 

where 
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s I 
C 


s f 
XK1  . 


t i 


stiffness  ol  tin-  axial  sprint 
stiffness  of  the  *;li* • i r sprint' 
damping  in  axial  load 
damn inc  in  shear  loading 
stiffness  of  restoring  torn ue  sprit 
damning  of  restorin',  torque  sarin 


And  the  subscript  "n'  denote^'  the  nuantitv  in  the  itpd'->'or'  od  'ate,  t 
t“0 , and  (')  denotes  differentiation  "Ith  respect  to  r-'me.  The  fc  e>  , 

T6Y  . , T6X.,  and  the  moment  B^+j  are  derived  hv  cl'aimln"  the  su!  scr!"t 
i on  the  right  side  of  equations  (2)-(4)  to  J+l. 

The  articular  facets  have  iieen  modelled  hv  two  sprints,  one  limit :r  ■ 
rotation  and  the  other  limiting  the  sliding,  of  one  vertebra  over  the 
ad-jacent  ones.  The  force  resisting  relative  rotation  het ue.  > the  ith  and 
the.  i-1  vertebrae  in  the  erect  mode  is  given  bv 

T51  = XKhL  • h^sin  f (0,^  - - (01  - 0^) } (5 

where  h.  = the  distance  of  tin  articular  f.  'w  * ’ e r o * ’ • 

i 

vertebral  bodv  and  XKh.  * the  stiffness  of  tie  snrtrg,  resistor"  the  m In- 
ti ve  rotation. 

Tn  t!ie  hyperextended  node,  due  to  the  forced  change  in  curvature 
of  the  spine,  it  is  assumed  that  the  facets  have  bottomed  out.  ar.  bene 
the  lamina  act  as  another  beam  parallel  to  the  dish.  The  c* ange  in  length 
AjA,;  as  shown  in  Figure  1°  is  used  as  t'>e  measure  of  axi-1  deformation  of 
this  beam,  and  the  force  developed  is  given  'v 

T51  - XKht  *(A1A2i  - A1A2f)()  < ■ 

The  force  resistinc  sliding  motion  at  t’'c  facets  "fau-  t 1 • ■ it1'  red 
the  i-1  vertebra  is  given  bv 

FX1-  VKhj  «(A2A31  - A2A10{) 

Due  to  the  overmanning  nat  ire  of  tl  e articular  facets  It  is  'ffJco 
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to  define  the  point  of  application  of  the  foregoing  force.  However,  in  the 
model,  it  lias  been  assumed  that  this  force  acts  perpendicular  to  the.  longitud- 
inal axis  of  the  vertebral  body  at  a point  (d.  + AC-j  /2)  below  the  center  of  the 
vertebral  bodv,  Figure  20.  The  reaction  of  this  force  on  the  vertebra  irmediat- 
lv  !e  low  acts  at  a distance  d.  , above  the  carter  of  the  i-1  vertebra.  This  as- 
sum; tion  is  justified  because  the  superior  articular  facets  are  shorter  than  the 
lamina  and  the  inferior  articular  facets.  As  shown  in  figure  20  the  forces 
on  the  itl  vertebra  parallel  to  un  and  w^-axis  are 


J1F  . = ('rhv,  + fv2  - FV1)  cose.  - (TbY.  + T52  - Til)  sinO. 

in  l i i i (8) 


- T7yt  cos  {e.  - (0oT0ol..1)  } + T7Y.  sin  (0.  - C - 0C  ±_1 ) } 

;:f  . = (T6Y  + T52  - T51)  cosfl . + (T6X . + FX2  - FXl)  sin0. 
ui  l ii  l 

- T7Y,  cos  (0.  - 001_j)  ) - T7X,  stn  (0,  - (<!„,  - 

e0i-l» 

Taking  the  sum  of  moments  about  the  center  of  gravity. 


;:v_.  = T7Y.  cos  fen.  - 0_.  ,}e.  + B.  , - T6Y.C.  - Tbx.d.  - T7Y . sin 
r,i  i 0i  Oi-l  i i+1  i i ii  i 

(e0i  - °0i~i>di  - T7yi  co*  {90i  - e0i-i}di  - T7xi  Bin  CIO) 

{e0l  - 0Oi-l}ei  - + CT51  - T52)(e.  + h^  - FX2.dJ  - fvl(d.  + AC.±/2) 

There  are  auxiliary  forces  and  moments  which  are  to  be  added  to  these 
equations  whenever  applicable.  Such  loads  Include  the  shoulder  strap  force, 

Ian  belt  force,  seat-back  reaction,  chin-chest  contact  forces,  and  reaction 
from  the  hyperextension  block.* 

Finally,  the  equations  of  motion  can  be  written  as 


7.  V , = 

IU 
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i 
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where  nij  is 

the  mass 

of  inertia  ; 

ilvmt  the 
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For  each  boJv  there  are  three,  such  equations,  hence  the  total  number  of 
equations  to  he  solved  simultaneously  for  i = 1—26  is  78.  This  "as  a.hiev  ’ 
bv  using  Hamming's  Predictor  Corrector  method.  T|)P  iccpi erat ' ons  ' and 
Y are  input  parameters  corresponding  to  t' e sled  ircelnratinn  in  v 

and  v directions. 

4.  Choice  of  Parameters. 

Ceonetrv  of  the  Spine.  The  model  requires  the  initial  c.oor di nates  r>r  t-',£ 
center  of  each  vertebral  bodv,  ti  e length,  the  angle  the  Ion  - i f u-’i  ml  rrw 
of  the  vertebra  makes  with  the  vertical,  the  distance  of  the  articular  facets 
from  the  center  of  the  vertebral  body  and  the  thickness  of  the  'isk  coincident 
witli  each  vertebral  bodv.  These  measurements  worn  made  tb  t1  e.  he1"  nr 
x-ravs  taken  before  each  run,  with  the  cadaver  placed  in  the  accelerator. 

Mass  and  Moment  of  Inertia  of  F.ach  pi  gid  I. ini  . ’the  model  requi  res  the 
distribution  of  the  mass  of  the  whole  body  to  the  rigid  .links,  kin  and 
Wickstrom  (34)  have  estimated  the  mass  distribution  of  segmented  cadaveric 
trunks.  The  mass  distribution  used  in  the  model  is  consistent  with  their 
data.  The  three  different  cadavers  used  were  weighed  before  the  experiment. 
After  the  experiments,  the  torso  was  cut  at  the  level  of  1.3  and  the  two 
halves  were  weighed  again  to  give  the  total  mass  above  1.3.  Three  further 
segments  were  made  of  the  unner  torso' — the  first  consisting  or  R3  to  T8, 
the  second  from  TP  to  C7,  and  the  third  r rom  r7  to  the  head,  the  f'.r-e  qponenf" 
were  weighed  again,  their  centers  of  gravity  - ere  -’otermined  sing  a.  Ira 
platform  which  utilizes  three  load  cells  at  fiat".)  distune  s fro;  i one  anoth"r, 

A trifilar  pendulum  was  used  to  determine  the  moment  of  inertia  of 
these  segments.  Based  on  the  foregoing,  measured  data  the  mass  and  moment 
of  inertia  at  the  various  vertebral  levels  were  estimated.  ’r!  e marts  of  r'  •> 
arms  was  assumed  to  be  distributed  to  the  first  f,-ve  thoracic  vertebrae. 

r“^~  " ■— 11  — 
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ri.vriic.il  Constants.  l’ach  rigid  link  in  the  model  is  associated 
ill  i disl  requiring  throe  .prii  g constants  aid  facets  requiring  two  spring 
<•  instants.  'once,  a total  of  131  spring  constants  are  required  for  the 
model.  Theoretically,  all  130  constants  can  be  different  and  in  an  ideal 
situation,  they  have  to  he  experimentally  determined  for  eacli  specimen. 

•rv  little  data  are  avialab le  in  the  literature,  and  in  the  case  of  the 
articular  facets,  nothin;;  is  available.  King  and  Vulcan  (35)  have  measured 
lie  stiffness  in  compression  of  the  lumbar  and  the  lower  thoracic  vertebrae, 
arkolf  (3b)  ’.as  measured  tie  axial  stiffness  of  disks  from  1,4  to  ^8  and 
as  found  values  ranging  from  7000  lb/ in.  in  the  lumbar  region  to  19000  lb /in. 

: n t'.  thoracic  region.  No  data  is  available  for  the  upner  thoracic  level. 

'arkolf  also  gives  the  rotational  stiffness  of  the  disks  from  700  in-lb/rad 
to  2400  in-lb/rad.  However,  these  measurements  were  made  with  no  preload 

- 

and  at  very  small  deflections.  Most  1 iological  materials  possess  a sig- 
moidal load-deflection  characteristic,  that  is,  the  stiffness  increases  with 
increasing  loads.  Hence,  the  foregoing  values  for  rotational  stiffness  of 
t . ■ disks  appear  to  be  much  too  low;  and  in  the  lumbar  region,  a rotational 
. tiffness  of  6000  in-lb/rad  was  assumed.  In  the  thoracic  region,  the  rib 
cage  imparts  more  rigidity  to  the  spine,  and  a rotational  stiffness  of 
12000  in-lb/iad  was  assumed.  In  the  cervical  region,  a rotational  stiffness 
of  2400  in-lb/rad  was  used  based  on  the  values  given  by  Vulcan  (37).  Typical 
constants  used  ir  the  model  are  shown  in  Table  7. 
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Table  7 List  of  Constants  Dsed  jn  Spinal  ’’odel  (Cadaver  2 li  i l ) 


Li  near 

r ot  itinnnl 

Vo t ;■  f ’ ft  ; 

Moment  of 

Damping 

c‘  t i f fness 

’'amu  ; n e 

Mass 

Cl  h-  sett  / 

1 ne  r I t a 

L \ nen  r 

St  if fn«- 

?s  Disc 

Of  sc 

0 i sc 

(lb- in- 

D I :r 

Facet 

(lb-sec/ 

On-P  / 

(in-!' 

Segment 

in) 

sec'  ) 

(lb/ in) 

(lb/ir ' 

in) 

rad ) 

sen /rn  '3 

Head 

0. 19722 

0.2816 

29991 

1200 

6.34 

2400 

? . 0 

ri 

0.01250 

0.001 3 

1595 

1200 

j . T«i 

3AOO 

1 . n 

C2 

0.01250 

O.0013 

1 595 

1200 

1.30 

2400 

1 .0 

C3 

0. 00441 

O.onift 

1595 

1200 

! . J'' 

2400 

1.0 

04 

0.00441 

0.0025 

1595 

1.200 

1 . 30 

2400 

1.0 

C5 

0.00441 

0.0030 

2537 

1200 

2.ftc 

2400 

1.0 

C6 

0.00441 

0.0030 

2537 

1200 

2.60 

2400 

1 .0 

C7 

0.00441 

0.0040 

2537 

1200 

2.60 

2400 

1 .0 

T1 

0.00441 

0.0066 

16562 

2000 

25.00 

12000 

2” . n 

T2 

0.01031 

0.0184 

16562 

2000 

25.00 

1 2000 

,30.0 

T 3 

0.01031 

0.0255 

16562 

2000 

25.00 

12000 

’>0,0 

T4 

0.01031 

0.0273 

26562 

2000 

25.00 

1 2000 

29.0 

T5 

0.01031 

0.0340 

16562 

2000 

25.00 

12000 

20.0 

Tf> 

0. 0165ft 

0.0392 

10310 

2000 

20.00 

12000 

20.0 

T7 

0.0165ft 

0.0476 

10310 

2000 

20.00 

1 2ono 

?o.n 

T8 

0.01656 

0.0491 

10310 

2000 

15.00 

12000 

30,0 

T9 

0.01656 

0.0546 

10310 

2000 

1 5 . n0 

1 2000 

20,0 

T10 

0.01656 

0.0537* 

8000 

2000 

1 5 . 00 

12000 

■>o  0 

Til 

0.00254 

0.0625 

8000 

2000 

10. 0 

6000 

30,0 

T12 

0.00254 

0.0156 

8000 

2000 

1 o , 00 

4000 

10,0 

LI 

0.00257. 

0.0156 

8000 

6000 

10.00 

6000 

10.0 

L2 

0.00160 

0.0156 

8000 

6000 

10.00 

( 000 

10.0 

13 

0.00160 

0.0156 

3000 

6000 

10.00 

6000 

10. 0 

L4 

0.00160 

0.0700 

8000 

6000 

10.00 

6000 

10.0 

L5 

0.00160 

0.0700 

8000 

6000 

10.00 

6000 

10.0 

Pelvis 

0.02991 

0.4540 

10000 

0 

240.O-) 

0 

20. 0 

Stiffness  and  damping  values  refer  to  the  disc  nr  facet  fcelo>;  tbe  '-erne.nl: 

5.  Model  Validation 

Fxneriments  involvin'’  the  use  of  human  cadavers  i.’ere  carried  out  to 
validate  this  model.  In  both  the  computer  an*'  exnerimertal  simii1  •’ tiorr  , ’•he 
conditions  at  ')Otb  ends  of  the  spinal  column  are  .specified.  'n  ac col er  ■ • inn 
fnout  is  applied  at  the  neJvfs  ■ iile  the  top  of  .bead  is  ’n  remain  trer< 
free  at  all  times,  ’’j nematic  a './or  kinetic  'uiantiti.es  can  1 measure  i t' 
vlipre  alonp  the  column  for  comr, orison  wi>  It  those  predicted  hv  the  model. 

The  parameter  used  r or  validation  in  tbi  ease  1:  the  force  1 too 


5" 

i ilia  vertebral  bodies.  An  ii  tervjrtebra)  m»Iw  load  coll  (T'/T.c)  was 
•va- :*>;>•  •’  tu  provide  (lie  maRnitiue  an  1 the  line  of  action  of  this  force. 
•:\!»eriii>  nfui  techniques  and  i ns  t.  rumen'  at  ion  used  h *ve  been  describe' 
in  detail  Hv  Prasad  et  a!  (18).  However,  it  should  he  pointed  out  that 
vie.  r 'i'  was  instrumental  in  the  establishment  of  th<  role  of  the  articular 
facet'  and  rfe  posterior  spina1  structures.  Facet  loads  were  obtained  bv 
'irst  computing  the  total  spine  load  at  the  level  of  the  IVLC.  Tt  was 
i . i ‘ f r ’-.e  the  product  of  t he  measured  seat  -pan  load  and  the  ratio  of 

i !•  • torso  weight  above  tie-  fVT.C  to  the  total  bodv  weight.  "'lie  difference 

■ ta  s'l.  tin  total  mine  ’o.ad  and  the  intervertebral  load  is  the  facet  load. 

■ • «.  installed  on  H'e  posterior  surface  of  the  lamina  provided 

» :-lit tonal  qualitative  evidence  of  the  facet  load.  As  shown  in  Figure  12 
, In*  fae.-t  load  curve  is  mirrored  by  that  of  the  posterior  strain  gage,  in- 

i udine  tlie  instant  of  time  at  which  a c bailee  in  sign  occured. 

■'dree  cadavers  'mre  subiected  to  -HI  or  seat-to--head  tranezoidal 
icceleration  pulses  of  6,  8,  and  10  q,  in  the  erect  and  hvperextended  mode. 

■ rate  of  onset  was  approximately  250  p/nec.  Figures  21-2^  show  compari- 
uns  >>r  each  of  the  three  cadavers  at  different  p levels  and  spinal  con- 
, orations.*  A comparison  of  both  the  intervertbral  load  and  the  facet 
t..r  ( a da  vc  r 2231  is  s'’own  in  Figure  ’1.  This  run  was  made  at  10  « in 
hi  e cct  me.:.;.  These  load:  were  measured  at  the  level  of  1.2-T.3.  llote 
the  change  in  sign  of  the  facet  load  as  a result  of  the  eccentricity  of 
tie  head  arid  torso  with  respect  to  the  spine.  In  both  the  model  output  and 

•; 1 1 experimental  lata,  the  secor.  1 intervertebral  force  peak  coincides  with 

» 

...  maximum  tension  in  the  facets,  Fi cure  22  shows  a run  in  the  hvnerextended 
iode  using  the  same  cadaver  and  the  same  level.  There  is  a marked 
* f'om:>arisnn  for  all  Is’  runs  are  p,ivon  hv  ^r.isnd  (in)  • 


_____ 
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reduction  in  intervertebral  load  arcorin-m i e*l  by  an  increase  in  compression 
at  the  facets  which  did  not  go  into  t'-nsion  or  unload  as  before.  'T’be  node1 
predicted  this  phenomenon  and  offers  an  explanation  for  the  dramatic  in 

crease  in  fracture  p level  reported  by  Fwinp  et  -a!  ( 4 ) . Tt  is  interest  fc" 

■ - 

to  note  that  there  was  an  increase  of  on]'-  110  1 in  the  peak  intervertebral 
load  when  the  g level  was  increased  from  8-10  g.  The  change  was  1156  11  for 
a similar  increase  from  6-3  g.  The  model  predicted  a chance,  of  334  " ror 
the  first  case  and  1121  N for  the  second. 

A 6 g run  on  Cadaver  2200  is  shown  in  Figure  23.  this  is  a connarison 
of  a run  made  in  the  erect  mode.  Again,  there  was  unloading  of  the  facet'-. 

A hvperextended  run  at  8 g on  Cadaver  2413  is  shown  in  Figure  24.  Al- 
though the  peak  values  of  the  intervertebral  loads  matched,  the  choice 
of  the  model  constants  resulted  in  a system  with  a higher  natural  frequence 
than  that  of  the  cadaver.  The  facets  barelv  went  into  tension  towards  the 
end  of  the  pulse  and  the  model  was  able  to  nredict  this. 

6.  Discussion  and  Conclusions 

A mathematical  model  of  the  spine  has  been  formulated  to  sirmlw,  ji 

response  to  Impact  acceleration  in  the  sagittal  plane.  h has  » ecu  vn'Mat'-,’ 

by  comparing  its  predictions  wit!  experimental  data  derived  fro-  1°  runs 

on  3 cadavers  which  were  suhiected  to  3 level;  of  -M'  arceler  U inn.  ''he 

z 

question  of  what  constitutes  proper  validation  of  a model  cannot  h ru!l  ■ 
answered.  However,  some  basic  ground  rules  wore  established  for  *-'m-  va ' - 
idation  procedure.  Since  the  acceleration,  levels  were  relatively  -v-d 

the  differences  were  small,  it  was  decided  that  there  vns  to  be  onlv  or< 


set  of  constants  for  each  cadaver  for  all  6 runs.  Variation  among  cadavers 
was  permitted  since  it  is  well  'Town  that  biological  differences  arc  cm'te 
large.  However,  the  values  of  these  constants  must  rot  deviate  too  far 
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from  t hose  published  in  fi  e literature.  for  example,  tht*  stiffness  of  the 
sprine  representing  the  lower  thoracic  spine  assumed  values  of  8000  to  10000 
lb/in.  for  the  3 cadavers.  Those  for  the  lumbar  spine  remained  at  8000  lb/in. 

In  fact,  the  same  linear  and  rotational  cl  amp Inf  coefficients  were  also  used 
for  all  3 cadavers.  In  my  case,  no  serious  attempt  was  made  to  optimize 
• constants  to  fit  tlu  experimental  data.  There  was  an  extreme  paucity 

| 

of  data  as  far  as  the  stiffness  of  the  facets  and  the  posterior  structures 
of  the  vertebra  is  concerned.  Although  the  model  provides  for  the  selection 
of  different  spring  rates  for  the  facets  in  compression  and  tension,  the  same 
value  was  used  as  a first  approximation.  It  can  be  argued  that  in  tension 
the  joint  capsule  and  all  posterior  ligaments  are  in  action,  whereas  in 
compression,  there  is  only  bony  deformation.  The  model  was  quite  stable 
in  terms  of  its  ability  to  accept  a wide  range  of  values  for  the  constants. 

The  numerical  computation  can  usually  he  carried  out  to  the  end  of  a rea- 
sonable; input  pulse  (200  ms)  without  any  breakdown  and  the  results  are 
almost  identical  with  different  time  steps. 

The  significance  of  the  initial  curvature  of  the  spine  was  brought  out 
bv  this  model,  as  evidenced  hy  the  difference  in  response  between  the  erect 
and  hyperextended  modes.  The  ability  of  the  model  to  predict  the  change 
in  sipn  of  the  facet  force  is  even  more  important  since  the  unique  feature 
of  a dual  load  path  is  the  key  element  in  a full  understanding  of  the 
mechanism  of  fracture  due  to  +0  acceleration.  The  anterior  aspect  of  the 

7 

vertebral  body  Is  subiected  to  an  extremely  high  compressive  load  when  the 
facets  unload  and  go  into  tensioi.  The  experimental  basis  of  this  phenomenon 
has  been  fully  documented  by  Prasad  et  al  (38),  and  this  model  is  proposed 
as  a candidate  for  an  inlurv  model  to  replace  the  present  one-degree-of- 
f reedom  base  excitation  model  which  uses  the  concept  of  dynamic  overshoot 


v 


>9 


f 
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as  a measure  of  peal,  load  on  the  spine*  in  fact,  there  is  vorv  little  dvn- 
amic  overshoot  in  the  cadaver  data  acquired  thus  fir.  "'I.  first  pea  in  fie 
intervertebral  load  is  invariably  lower  than  the  - erond  ror  a run  made  in  ’lie 
erect  mode.  That  is,  the  noal  load  e;:ner  owed  hv  tin  vert  «•>  i 1 bodies  ncci" 
much  later  in  the  pulse  and  is  a result  o!  head  apt  n.r*„  ! inenatics  ,pd 
load  transfer.  Tt  is  not  the  result  of  a dynamic  overshoot. 

VIII.  1TECHANISM  OF  TN.TttP.Y  AND  INJURY  PREDICTION 

On  the  basis  of  the  data  obtained,  it  is  now  possible  to  assemble  the 
research  results  on  spinal  iniurv  and  nronosc  an  iniurv  "-echanirn  that 
fully  explains  the  common 1\  observed  anterior  wedge  fracture-  sustained  W 
pilots  durinr  emergency  egress.  Previous  worl  hv  ’’ule.i-  nd  nr  (39'  es- 

tablished the  fact  that  during  caudocenhalad  (4-r  ) nccel  °r«t -‘i  a , the  spine 
is  subjected  to  both  axial  compression  and  bending.  The  bendity  effects 
are  due  to  the  eccentricity  of  the  torso  with  resnect  to  the  snine  and 
are  enhanced  bv  the  forward  rotation  of  the  head  and  torso.  Subsequent 
work  by  Ewing  et  al  (A)  showed  that  fracture  levels  could  bp  raised  sig- 
nificantly bv  placing  a 5.7  cm  thick  hvperextension  block  opposite  I.J  on 
the  se.nt  back.  The  curvature  of  the  spine  was  altered  !"•  the  1 lock,  but 
there  was  insufficient  movement  to  decrease  the  eccentric! t'  of  the  torso  and. 
hence,  the  bending  moment  on  the  spine.  Nevertheless,  tr’-o  average  level 
of  fracture  wn«  raised  from  1 0 tc  18  g,  and  M - de  rease  ;n  . interior  strain-’ 
was  generally  significant  at  the  *5  percent  con1  idem-  level.  rbe  b"PO-- 
thesis  that  the  facets  act  as  not  ion  limiters  '<-d  to  this  docunontat  ion  of 
their  role  during  +<  ^ acceleration.  In  the  erect  spiral  mode  t e facets 
tend  to  unload  and  go  into  tension,  causing  the  vertebral  bodies  to  aufta'n 
more  comnressive  load  than  the  total  snine  load,  "ids  occurs  wV-n  tie  bead 


an-!  tors  undergo  maximal  fer-  ir'  flexion.  T1  --  u • -r’n  -r  •-  , 
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,ro  therefore  the  result  of  eccentric  compression  coupled  with  the  unload- 
ing of  t!u:  facets.  T.n  the  hyperextended  node  the  facets  relieve  the  verte- 
’ r-tl  bodies  of  some  of  t o compr'  ssive  load,  and  thus  i t was  nossible  to 
raise  tie  fracture  level  by  such  a considerable  margin. 

The  model  output  cue  It  used  to  predict  the  forces  and  moments  sus- 
tained by  each  vertebra  at  each  instant  of  tine.  The  intervertebral  axial 
h.irce  is  plotted  as  a function  of  distance  above  the  seat  pan  in  Figure  25 
for  a in  r run.  mhe  curves  are  plotted  for  instants  of  time  at  which  the 
first  and  second  peaks  occur  and  for  both  the  erect  and  hvperextended  mode. 

"ote  that  the  intervertebral  load  is  less  in  the  hyperextended  case  than 
that  for  the  erect  mode.  Also  the  lower  thoracic  spine  sustains  a pro- 
portionally higher  load  and  is  therefore  more  susceptible  to  fracture.  These 
"uns  simulate  the  cadaver  runs  carried  out  on  the  vertical  accelerator,  and 
since  the  cadavers  were  eviscerated  for  the  installation  of  strain  pages 
and  the  intervertebral  load  cell,  the  inertial  load  on  the  lumbar  spine  was 
minimal . T>v  use  of  the  mass  distribution  data  provided  by  Liu  and  Viekstrom 
(34),  it  was  found  that  the  lumbar  spine  sustained  a slightly  higher  load 
at  both  the  first  and  second  peaks. 

If  the  dvnamfc  fracture  load  for  vertebrae  is  known,  this  model  can  be  used 
is  an  injurv  predictor.  The  fracture  loads  can  be  plotted  as  a band  on  figure  25 
,u.J  inlur”  is  1 ikelv  if  the  force  sustained  bv  anv  vertebra  exceeds  its  fracture 
level  as  a conservative  estimate,  the  static  fracture  loads  given  bv  Ruff  (7), 
can  be  used  as  a limit,  pending  the  acquisition  of  dynamic  data. 

IX.  CONCLUSIONS 

Sufficient  data  are  now  available  to  formulate  a comprehensive  mechanism 
O Injury  to  the  snine.  A mathematical  model  has  been  developed  and  vali- 
dated. It  can  be  used  to  further  the  understanding  of  this  mechanism  and 
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ro  predict  injurv  as  a result:  of  an  arbitrary  input  pulse.  Since  the 
mechanism  is  not  one  of  simple  c impressive  failure,  it  was  found  that 
t t restraint  .>  item  presently  available  in  military  aircraft  can  l e used 
feet  ive.i  v to  h nerextencl  the  spine  and  thus  raise  the  fracture  g level 
•v  a considerable  margin.  he  existence  <f  a dual  load  oath  alone  tin 
se'iic,  one  ••  .i>  ««<»'•  the  vertebral  bodies  and  the  other  along  the  facets , and 
posterior  structures,  has  been  documented.  As  a result,  a logical  ex- 
of  -nation  < in  he  given  for  the  commonly  observed  anterior  wedge  fractures. 

''he  data  acquired  from  canine  experiments  were  incorporated  into  a 
•iodvnanic  model  of  the  spine  which  was  developed  under  another  OMR  Contract 
o.  1000 14-75-C- 1011.  Preliminary  results  of  this  model  were  given  hv 
Tennyson  and  King  (40).  T ts  output  was  validated  against  human  volunteer 
data  acquired  by  Twine  and  Thomas  (41). 

!The  results  of  this  research  program  can  be  used  hv  the  Haw  to  improve 

the  design  of  ejection  seats,  existing  seats  can  also  be  modified  to 
it. elude  spinal  i vperextension  as  part  of  the  ejection  sequence.  Present 
technology  is  capable  of  effecting  hyperextension  in  a short  period  of 
time  prior  to  ejection. 

future  research  should  be  directed  towards  the  development  of  a sur- 
rogate sp’ne  which  possesses  human-like  responses  to  +C^  acceleration. 
'inenatic  data  of  spinal  motion  are  required  to  define  corridors  of  spinal 
response  to  axial  compression  as  well  as  flexion  and  extension.  High 
precis  ton  techniques  in  accelenmetry  will  be.  required  to  achieve  this  end, 
'Ain  use  of  such  surrogates  along  with  mathematical  models  constitutes  an 

« 


ideal  approach  to  the  design  of  safe  ejection  seats  of  the  future. 
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